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Dynamical Fermion on Random-Block Lattice 
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Massless fermion field interacting with abelian dynamical gauge field on 2-dimensional random-block lattices 
are investigated using Hybrid Monte Carlo simulations. Preliminary results of the Wilson loop and the chiral 
correlation function are in agreement with the continuum Schwinger model. 



1. INTRODUCTION 

The problem of formulating the chiral fermion 
field on space-time lattices has been discussed 
from the viewpoint of funtional integrals Q. I 
argued that the functional integral measure pro- 
vided by any single lattice ( hypercubical lat- 
tice or random lattice ) is not sufficent to yield 
fcrmionic functional integrals to agree with the 
continuum field theory. My proposal of using 
an ensemble of random-block lattices ( RBL ) to 
provide a proper measure has been investigated. 
RBL regularization has been tested for the free 
fermion ( massive and massless ) in 2d and 4c? [^[U 
respectively, and for 2d massless fermion in an ex- 
ternal abelian gauge field [[|J|] . All observables in- 
cluding fermion propagators, gauge invariant cur- 
rents, axial anomaly and current-current correla- 
tions are all in good agreement with the contin- 
uum field theory. It is natural for us to proceed to 
investigations involving dynamical fermion field 
on RBL. In this paper, massless fermion field in- 
teracting with the dynamical abelian gauge field 
on 2d RBL is studied. Using Hybrid Monte Carlo 
simulations, the Wilson loop and the chiral corre- 
lation function are measured and preliminary re- 
sults of these observables are in agreement with 
the continuum Schwinger model. More exten- 
sive measurements involving higher statistics and 
larger lattices are now in progress and will be re- 
ported elsewhere. 



2. SCHWINGER MODEL on RBL 

The action of a massless fermion field interact- 
ing with the dynamical abelian gauge field on a 
2d RBL is 

Arbi = WjK^j^fj, [U/j, {i)ipi+fi - U^{i - fytpi-p,] 



+(3j2w p [l-Re(U p ) 



where ipi and 'ipi are two independent 2- 
component spinors at site i, u>i — FT fe±±£zfEi=flj 
is the site weight at site i, = (xi + fc — 
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is the inverse of the distance between the sites 
i + /} and i — (x and U^(i) = exp[ie J x, +' 1 A ■ dl] 
is the link variable pointing from the site i to the 
site i + /t, and the 7 matrices are chosen to be 
7i = 01, 72 = ^2 and 75 = -«7i7 2 = cr 3 . The 
last term in the action A^i is due to the pure 
gauge field and is conveniently labeled as 

A g [U] = ^v3 p [l-Be(U p )] 
p 

where j3 = Jr, U p is the path-ordered product of 
the link variables around a plaquette p, w p = -g- 
is the weight of the plaquette p and S p the area of 
the plaquette. The action A r w is invariant under 
the Uy(l) transformation 

ipi — > exp(— iedijipi 

-tpi — > ipi exp(ie6' i ) 
U^i) -> exp(-ie0 i )£/ M (i)exp(ie0 i+A ) 
and the Ua{^) transformation 

ipi -> exp(ie7 5 a)^i 
ipi -> ipi exp(ie7 5 a) 
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where a is a global parameter. In the continuum 
limit, A r u goes to the Schwinger model || 

As = J d 2 x feb,4^ + ieA M (x)]V(a;) + 
The partition function of A r u is 
Z = /[dC/][#][#]exp(-A rW ) 
[d[/]det(G)exp(-A g [C/]) 
where G is the matrix 

Off = ifuiKRUMSjM Ul(i £)<^_ A ] 

and det(G) is positive definite even for one single 
flavor of the fermion field . With 71 = oi, 72 = 
02, G can be written as 



G 





Ti + iT 2 



Ti - iT 2 




where T\ and T 2 are anti-hermitian matrices, 
= -Ti and T 2 f = -T 2 . Then 



G 



and 



_(Ti + ir 2 )t 

Ti + iT 2 



det G = det(Ti + *T 2 ) • det(T x + iTa)^-!)" 4 
= |det(T! + iT 2 )| 2 
= det^if) > 

where _ftT = T\ + iT 2 and n s is the total num- 
ber of sites which is assumed to be even. There- 
fore the probability distribution function is posi- 
tive definite for any gauge configurations and thus 
amenable for Monte Carlo simulations. Complex 
scalar fields {<fii,<fi*} are introduced to convert 
det(G) to a functional integral over the complex 
scalar fields. The partition function becomes 

Z - / [««7][d0t][d0] exp (-A g [U] - cj)\K^K)- l 4,) 
\dU][^M]e X p{-A g [U]-CH) 



where £ = (K<) 1 <fi. Then the system is equiva- 
lent to that of effective action A e ff = A g [U]+^£ 



in which quantum expectation value of any ob- 
servables can be computed by the Hybrid Monte 
Carlo ( HMC ) simulation JtJ. The procedures 
of performing the HMC simulation are similar to 
those on a hypercubical lattice. 

3. OBSERVABLES 

Measurements of the Wilson loop and the chiral 
correlation function at (3 = 10.0 were performed 
using an ensemble of 64 random-block lattices 
each of size 12 x 12 and average lattice spacing 
a = 1.0. On each RBL, 2000 sweeps are used 
for thcrmalization, 1000 measurements were per- 
formed with consecutive measurements separated 
by 10 sweeps. More extensive measurements with 
higher statistics and on larger lattices are now in 
progress and will be presented elsewhere. 

CHIRAL CORRELATION FUNCTION 

One remarkable feature of the continuum 
Schwinger model is that all gauge-invariant quan- 
tities can be computed exactly. The most inter- 
esting physics of the model is that the chiral sym- 
metry is dynamically broken by quantum correc- 
tions. This can be observed in the chiral correla- 
tion function 



C(x,y) =< i/>(x) 



1 



r^-^{x)^{y)— 



•75 



which goes to a non-vanishing constant for large 
\x — y\, directly indicating the breakdown of clus- 
tering and degenerate vacua in the theory. In a 
finite square of size L x L, C(x, y) can be evalu- 
ated in the following. 



C(x,y) = 



exp[47r(A (a;, a:) - A(x, x) - A (x, y) + A{x, y))} 

where S fl (x 1 y) are components of the mass- 
less fermion propagator, Ao(x, y) is the mass- 
less scalar propagator and A(x,y) is the massive 

scalar propagator with mass \ f^- They can be 
expressed in the following. 

S F (x,y) = nfpSpfay) = 7 M 9 M Air(a;,y) 
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where 



1 

L2 



E 



(2ni + l)7r i 2 
ni,n2=— oo L L J 



|- (2n 2 + l)7r j 2 



ni ,n2— — oo 



exp[i ^-n- ( x - y)] 



A (^y) = i2 E 



ni ,ri2 — — oo 



2_7l27T \2 

exp[i^-n • — y)] 



On a RBL, the chiral correlation function can be 
measured as 
1 



Q{x,y) = -j [dU] det{G) e X p(-A g [U])x 

= ±j [dUm][dg]\K-\x,y)\*eM-A g [U] -C f 

where K~ l {x,y) is the right-handed fermion 
propagator in background gauge field and com- 
plex scalar field, which can be computed in the 
following. 



where Q = K^K and the last step is computed 
by the conjugate gradient method. The ensemble 
average is obtained by repeating the same calcu- 
lation for each RBL while holding the end-points 
x and y fixed in all RBL, 



< C(x,y) > 



n 



8(z 1 - x)S(z 2 - y)C(x 7 y) 



N, 



I alt 



E 



i=i 



This RBL regularized chiral correlation function 
< C'(x, y) > is then compared to the continuum 
result. In Fig. H, < C(x, y) > is plotted for an 
ensemble of 64 RBL each of size 12 x 12 and av- 
erage lattice spacing a — 1.0. It agrees with the 
continuum Schwinger model and goes to a non- 
vanishing constant at large \x — y\. 




Figure 1. The chiral correlation function 

< C (x,y) > plotted vs. \x-y\ at = 10.0. The 
solid line corresponds to the continuum Schwinger 
model. 



WILSON LOOP 

The chiral symmetry breaking in the Schwinger 
model results in complete screening of the lin- 
ear confinement potential of external charges 
in the pure gauge theory, and the quan- 
tum expectation of the Wilson loop W(C) = 

exp (^ie § c A ■ dx^j would display asymptotically 

perimeter law rather than the area law. In con- 
tinuum, 



< W(C) >= exp 



~4^J f d x ndy' i Ko(m\x-y\) 



where m = y —■ On a RBL, the Wilson loop can 
be measured as 



W(C)i = 



J[dU}[d<p}[diP}W(C) exp(-A rbl ) 
/[dff] [#][#] exp(-A r&i ) 
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_ j[du]m[de}w(c)exp(-A g [u] - eo 

J[dU][dt;}{^]exp(-A g [U}-eO 

The ensemble average can be obtained by repeat- 
ing the same calculation for each RBL while hold- 
ing the boundary of the Wilson loop fixed in all 
RBL. For a square Wilson loop of size D x D, in 
practice, we can increase the statistics by mea- 
suring all Wilson loops having D links in both 
directions respectively, and the average area of 
these Wilson loops is D x D. Unlike the chiral 
correlation function or other observables involv- 
ing fermion fields, the fluctuations of the Wilson 
loop from one RBL to another is rather small and 
the number of RBL used for ensemble averaging 
can be minimum. In Fig. ^, the quantum expecta- 
tion values of D x D Wilson loops for an ensemble 
of 64 RBL are plotted vs. D. The values are in 
good agreement with the continuum Schwinger 
model, indicating the complete screening of the 
confinment potential in the pure gauge theory. 

4. CONCLUSIONS and DISCUSSIONS 

The present investigation of the dynamical 
fermion field on RBL is far from completed yet. 
Besides the Wilson loop and the chiral correla- 
tion function, other observables must also be mea- 
sured, all must be with higher statistics and on 
larger lattices as well as for many different (3 val- 
ues. These measurements are now in progress 
and the results will be presented elsewhere. The 
present preliminary results are consistent with 
the conclusions drawn from previous studies of 
free fermion and fermion in a background gauge 
field on RBL, that an ensemble of RBL could pro- 
vide a proper measure for functional integrals in- 
volving fermion fields. Without breaking the chi- 
ral symmetry at the tree level, RBL regularization 
uses the naive fermion which presumably suffers 
from the species doubling on any one of the lat- 
tices, however miraculously gives the correct re- 
sult when the functional integrals are sumed over 
all RBL. The most remarkable feature is that 
it produces the correct chiral symmetry break- 
ing through quantum corrections. Whether RBL 
regularization also works for the chiral guage the- 
ories, say chiral Schwinger model or the Elec- 
troweak theory is still an open but very inter- 
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Figure 2. The values for D x D Wilson loop 
plotted vs. D at (3 = 10.0. The solid line corre- 
sponds to the continuum Schwinger model while 
the dotted line to the pure gauge theory. 



esting question. 
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